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The present article discusses magnetic confinement of the Dirac excitations in graphene in presence 
of inhomogeneous magnetic fields. In the first case a magnetic field directed along the z axis whose 
magnitude is proportional to 1/r is chosen. In the next case we choose a more realistic magnetic 
field which does not blow up at the origin and gradually fades away from the origin. The magnetic 
fields chosen do not have any finite/infinite discontinuity for finite values of the radial coordinate. 
The novelty of the two magnetic fields is related to the equations which are used to find the excited 
spectra of the excitations. It turns out that the bound state solutions of the two-dimensional 
hydrogen atom problem are related to the spectra of graphene excitations in presence of the 1/r 
(inverse-radial) magnetic field. For the other magnetic field profile one can use the knowledge of 
the bound state spectrum of a two-dimensional cut-off Coulomb potential to dictate the excitation 
spectra of graphene. The spectrum of the graphene excitations in presence of the inverse-radial 
magnetic field can be exactly solved while the other case cannot be. In the later case we give the 
localized solutions of the zero-energy states in graphene. 
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I. INTRODUCTION 

In recent years graphene, which is a single layer of car- 
bon atoms on a honeycomb lattice, has generated a lot of 
interest in the research community [lH3- This is essen- 
tially because of the experimental realization of graphene 
sheets in various laboratories and its possible applications 
in nanoscale devices. An electron in a graphene layer is 
described by a massless two-dimensional (2D) Dirac-like 
equation which produces a gapless linear energy spec- 
trum close to the Dirac points K and K'. The Dirac elec- 
tron imitates various relativistic phenomena like Klein 
tunneling 0] , unconventional Hall effect 0, d, [j| and Zit- 
terbewegung Q. To make electronic devices, the first 
thing is to confine the electrons. This can be attempted 
by either using an electric field or a magnetic field. Al- 
though zero-energy bound states or quasi-bound states 
may be created using electrostatic barriers or external 
potentials, the Dirac electron can not be always confined 
in this way due to its massless nature [7j. Particularly, 
the electrostatic barrier becomes completely transparent 
for normal incidence of the electrons. 

Some schemes have been proposed to tackle the prob- 
lem of confining the electrons. For example, A. de Mar- 
tino et al., has given an interesting proposal to confine 
the electron based on an inhomogeneous magnetic field 
Q. Later on various inhomogeneous magnetic field con- 
figurations have been used in an attempt to confine the 
electrons or create bound states @. Among the vari- 
ous field configurations, exponentially decaying magnetic 
fields [l0| , Gaussian magnetic field [ll[ , square well mag- 
netic barrier Q, non-zero magnetic fields in a circular 
dot [lH, magnetic fields corresponding to various solv- 
able potentials [lU etc., have been considered. In all 
these studies, discontinuous and/or inhomogeneous (in 
some specific direction) magnetic field profiles were con- 



sidered. 

Recently, Masir et al. Q has considered a finite-size 
magnetic dot with a constant magnetic field (B z ) within 
a circular region with radius i?, otherwise zero i. e. 
B = B Z Q(R — r)z. It is found that bound state does 
not exist in this configuration. In the present article 
our objective is to use different types of axial magnetic 
fields, which are continuous for any finite radial distance, 
to examine the possibility of creating bound states in 
graphene. In the first case we will consider the inverse- 
radial field pointed along the negative z-axis and whose 
magnitude is ~ 1/r. One can exactly solve for the Dirac 
electron-like excitations in graphene in presence of such 
a magnetic field. We will focus on the localized exci- 
tations in the present case. One of the surprising out- 
comes of the above analysis is an interesting relation 
between the bound state spectrum of the 2D Hydrogen 
atom and the energy spectrum of the magnetically con- 
fined graphene excitations. In this article it will be seen 
that the Schrodinger equation of the 2D Hydrogen atom 
plays a very important role for determining one of the 
components of the Dirac-like excitations of graphene in 
presence of the inverse-radial magnetic field. Due to the 
complex nature of the Dirac equation it turns out that 
the relevant Coulomb potential in the 2D Hydrogen atom 
problem becomes angular momentum dependent. Except 
this complication of an angular momentum dependent 
Coulomb potential the 2D Hydrogen atom bound state 
spectra gives us a clue of the bound state energies of the 
graphene excitations in an inverse-radial magnetic field. 

Although for an inverse-radial magnetic field one can 
find the spectrum of the localized massless excitations 
of graphene exactly, still the magnetic field profile is not 
realistic near the origin. Noting this fact, in the second 
case, we try to see the properties of the localized exci- 
tations in graphene in a more realistic magnetic field by 
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introducing a cut-off parameter (i.e, the magnetic fieid 
does not blow up at the origin). As a consequence of 
introducing this cut-off parameter the problem does not 
remain an exactly solvable one. Nevertheless localized 
zero-energy states can be exactly determined and it will 
be seen that in this case also the zero-energy states are 
infinitely degenerate. 

Before we list the materials of the present article in 
the last paragraph it is useful to give a brief discussion 
about the magnetic field profiles used in the article. The 
authors are aware of the fact that the kinds of magnetic 
fields used for the calculations are not generally used in 
laboratories. But in the future one may produce a mag- 
netic field which has a high but finite strength in the 
axial direction and slowly the field strength falls as one 
moves outwards from the axis. Using various cylindrical 
dipolar magnets, with the magnetic axis along the cylin- 
drical axis, with varying magnetic strength one can in 
principle produce a field which can resemble the realistic 
magnetic field with a cylindrical symmetry as assumed 
in this article. The central solid cylindrical magnet with 
the highest magnetic strength must be concentric with 
many annular cylindrical magnets with lower and lower 
magnetic strength. The magnetic strengths can be im- 
parted suitably in the laboratories. With two such set 
ups placed one above the other, such that the opposite 
poles face each other one can produce magnetic fields of 
various cylindrical symmetries in the intervening region. 

The materials in this article are presented in the fol- 
lowing manner. A brief and concise scheme of the mathe- 
matical framework which will be utilized throughout the 
article is presented in the next section. The spectrum of 
the localized excitations of graphene in the presence of 
an inverse-radial magnetic field is given in section UTTl In 
section IPV1 we present the exact zero-energy solutions for 
the graphene excitations in presence of a more realistic 
magnetic field which does not diverge at the origin. The 
last section IVl concludes the article with a brief discussion 
and summarizes the findings presented in the article. 



II. THE MODEL 

The quasi-particles in graphene follow a Dirac-like 
Hamiltonian given by 



H = V F (T ■ P . 



(1) 



where vp w 10 6 m/s is the Fermi velocity, er = (a Xl a y ) 
are the Pauli matrices and P = p + A is the canoni- 
cal momentum with vector potential A. For the sake of 
simplicity as well as clarity we will use the convention in 
which h = e = 1, where e stands for the magnitude of the 
electronic charge. The eigenvalue equation for the above 
Hamiltonian is Hip(x,y) = Eijj(x,y) where ip(x,y) is a 
two-component spinor given by 



■02 (x, y) 



Here the two components ipi^ix, y) denote the electron 
amplitude on two sites in the unit cell of a honeycomb 
lattice. We now write the eigenvalue equation in the 
following coupled form as 



n_^2 = eipx , n + V>i = eip 2 



(2) 



where e — E/vf and H± = P x ± iP y . The general form 
of the vector potential is of the form 



A, 



yf(r) ■ 



(3) 



where the function f(r) will be specified later. For this 
choice of the vector potential the associated magnetic 
field is perpendicular to the plane and it is given by 



B, 



-2/(r) - r/'(r) . 



(4) 



Since the magnitude of the magnetic field only depends 
on r it is convenient to use the plane polar coordinates 
in the present case. In plane polar coordinates one can 
write the component spinors as 



fa =r- 1 ' 2 e lme Xl {r). 



fa = r- 1 / 2 e^ m+1 ^ X2 (r) (5) 



where m — 0, ±1,±2, • • •, stands for the angular mo- 
mentum quantum number of the pseudo-spinor solution 
of the Dirac equation. Using the above solutions and 
Eq. ^ the eigenvalue equation for \i can be written as 



+ Vf(r) 



dr 2 



where the effective potential term is 
Vf{r) =Vi(r) 



1 1 



and 



V 1 (r)=r 2 f 2 -2mf-2f-r.f, 



(6) 



(7) 



(8) 



where the other term in the potential (r) is due to 
the centrifugal barrier. The factor of l/4r 2 appears be- 
cause of the presence of r -1 / 2 in the solutions in Eq. ([5]). 
Similarly for \2 we obtain 



dr* 



+ Vr(r) 



e 2 X2 ■ 



where the effective potential is 

(to + \){m + 



Vf{r) = V 2 {r) + 



and 



V 2 (r) =r 2 f 2 -2mf + rf. 



(9) 



(10) 



(11) 



It is easy to see that the above equations are standard 
radial Schrodinger equations (in 2D) with the effective 
potentials V[ c %{r). 
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III. THE INVERSE-RADIAL MAGNETIC FIELD 
AND THE SPECTRUM OF THE EXCITATIONS 

In this section we will concentrate on a magnetic field 
profile which is given by 



B- 



(12) 



where A > is a constant with the dimension of inverse 
length. In normal units the magnetic field will be rep- 
resented as —hX/er. A field of the above type can be 
obtained by choosing the /(r) as 



m = - . (is) 

r 

With the above form of f(r) one can see that 

(14) 



It is interesting to note that with the form of f(r) as 




FIG. 2: Plots of the effective potential V 2 vs y for m = 
(dot-dashed), m = 1 (dashed), m = 2 (solid) and m = 3 
(short-dashed). 



^ -1 




present case is related to the fact that the potentials de- 
pends on the angular momentum of the excitations. This 
observation is pivotal for the analysis which follows. 



A. The 2D Hydrogen atom and the upper 
component of the graphene pseudospinor 

To find the ground state and excited state energies and 
the corresponding wave functions of the Dirac electron- 
like excitations in graphene in presence of an inhomo- 
geneous magnetic field, as given in Eq. (|12j) . we require 



FIG. 1: Plots of the effective potential Vi S vs y for m = 
(dot-dashed), m = 1 (dashed), m — 2 (solid) and m — 3 
(short-dashed). 



given in Eq. (fT3|) . both the potentials Vi 2(f) turns out 
to be of the Coulomb- type for m > 0. The shape of 
the effective potentials Vi%(r) for some positive values 
of m are plotted in Fig. [T] and Fig. [5J The shape of 
the effective potentials clearly shows that there can be 
bound states for positive values of the angular momen- 
tum quantum numbers. For m < the potentials Vi^{r) 
are not attractive and hence do not support any bound 
state. Consequently, for states with m > the Dirac 
electron-like excitations in graphene due to the inverse- 
radial magnetic field can be understood by solving the 
radial Schrodinger equation, Eq. ©, in presence of a 2D 
Coulomb potential. The novelty of the 2D Coulomb po- 
tential, appearing in the Schrodinger equations, in the 



to solve Eq. 
given in Eq. 
^function 



where V\{r) in the effective potential is 
T4l) . One can write Eq. (|6]) in terms of the 



Mr) 



-1/2 



Xi{r) , 



dy 2 



1 d(j>x 
V dy 



2m + 1 

y 



y 



(15) 



0i=^i, (16) 



where 

y = \r, 

is a dimensionless parameter and 
A 2 £ = e 2 - A 2 . 
The energy eigenvalues are given by 



E = ±\v F \flTl- 



(17) 



(18) 



Eq. (|16|) can be converted to a confluent hypergeometric 
equation in x where 



V 



1 



P. 



N 



2m + 1 

N 
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for a real number N. Defining a new function <f>i(x) as 
fc(a:) = x m e- x /H x (x) , 

Eq. p^|) transforms into the following confluent hyperge- 
ometric differential equation: 



dx 2 



[(2m + l)-d^fi-(-iV- 



i)*x=o. 



(19) 



The most general solution of the above equation is 

$l(x) = A 1 F 1 (—N + rn + i, 2m + 1, x) 

+ B U(-N + m + ^, 2m + l,x), (20) 

where iF^-N + m + |,2m + l,a;) and f/(-JV + to + 
i, 2to+1, x) forms a pair of linearly independent solutions 
of the confluent hypergeometric equation. A and J3 are 
two arbitrary constants. 

In our analysis of 4>\{x) the angular quantum number 
to can be zero or positive integers. The nature of N will 
dictate the properties of the confluent hypergeometric 
functions. Out of various possibilities of choosing N one 
is particularly interesting. If N arc half-integers as 



2 ' 2 ' 2 ' 



(21) 



then — N+m+^ is a negative integer or zero if we restrict 
the to values in an appropriate way. Defining 



n = W+i = l,2,3,--., 

the possible to values for which — N + m 
negative or zero turns out to be 



(22) 



\ can be 



TO 



0,1,2,- 



1 



(23) 



When — N + m + h or 



+ to + 1 is a negative in- 
teger or zero then iF^— n + m + 1, 2to + 1, x) reduces 
to a finite polynomial in the positive x range. On the 
other hand if — n + to + 1 is a negative integer then 
U(— n + m + 1, 2to + 1, x) and iFi(— n + jn + 1, 2m + 1, x) 
are not linearly independent of each other and conse- 
quently the pair of functions in Eq. (|20p does not give 
the most general solution of Eq. (|19p . In this special 
case one may think of the second solution of Eq. (fT9|) as 
e x U(m + n, 2m + 1, —x) [l7| . which is linearly indepen- 
dent of iF 1 (— n + to + 1, 2to + 1, x). But, the function 
e x U(m + n, 2m + 1, — x) diverges at large values of x and 
consequently this second solution is also unsuitable in the 
most general solution of Eq. (pT9]) . 

If the second solution of the confluent hypergeometric 
equation is omitted then the solution of Eq. (|19l) coincides 
with the solution of the 2D hydrogen atom [18| . In this 
case 

x [ 1 F l {-n + m + l,2m + l,p n y)],{2A) 



up to a normalization constant and 
2(2to+ 1) 



2n- 1 



(25) 



The eigenvalues in the present case given by £ n>m will be 
of the form 



(2m + l) 2 
(2n- l) 2 ' 



(26) 



where m = 0, 1, 2, 3, • • • and n = 1, 2, 3, • • • the positive 
integers only. From the expression of the energy eigen- 
values as given in Eq. (|T5]) we get 



En 



(2to + l) 2 



(2n - l) z 



(27) 



The negative energies in the above represent the holes in 
the valence band. The exact energy eigenvalues of the 
Dirac electron-like excitations in graphene in presence of 
an inhomogcncous magnetic field is seen to depend on 
two quantum numbers. One of them is the radial one, 
n, while the other is related to the angular momentum 
quantum number to. 




FIG. 3: Plots of the normalized probability density of the 
zero-energy states vs y for n = l,m = (dot-dashed), n = 
2, m = 1 (dashed), and n — 3, m = 2 (solid). 



B. Degeneracies 

If one looks at the energy spectrum of the excitations 
one will see that there are states which are degenerate. 
The degeneracies can be understood from the form of the 
energy eigenvalues as given in Eq. ()18[) . It can be easily 
observed that 



Ei.o — E 2 .i — E-, 



3,2 



En 



0. 



(28) 
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This implies that the ground states of the excitations are 
infinitely degenerate. The corresponding wave function 
is given by 



V> (m+1 ' m) (r,6>) = N„ 



m e -y e imb 



(29) 



where m = 0, 1, 2, • • •. The normalization constant is 



N m = 



27r2 2 ( m + 1 )r(2r 



(30) 



The probability density plots for some of the zero-energy 
solutions are given in Fig. [3J All of the curves show that 
the states are localized near the origin and can be treated 
as proper bound states. 

In general, there are many other infinitely degenerate 
states. For example, the energy levels are infinitely de- 
generate if (2m + l)/(2n — 1) = p/q, where p and q are 
odd integers with p < q. The infinitely degenerate zero- 
energy state is obtained for q = p. 



C. The lower component of the pseudospinor 

For excited states, e ^ 0, the lower component of 
the pseudospinor can be obtained by the relation ip2 — 
in + i/>i. In the present case ipi(y) — e" ne (f)i(y) and con- 
sequently one can write 



e 



d to 

1 + 1 

ay y 



My), (3i) 



where <fii(y) is as given in Eq. (pM|) . But the lower com- 
ponent of the zero-energy states cannot be obtained from 
the formula given above. To obtain the the lower com- 
ponent of the pseudospinor for the zero-energy case one 
can directly use the first order equations in Eq. Solv- 
ing the first order equation Tl-ip2 — where II = 
e-' ie [-id/dr~(l/r)d/de + i\} and fa{r) = 2 (r)e l < m+1 ) e 
one easily obtains 



<fo{y) = 



y 



m+l 



(32) 



D. Reversing the direction of the magnetic field 

It is interesting to see what happens to the spectrum 
of the excited states when the magnetic field direction is 
reversed. From Eq. (TT31 it can be seen that the magnetic 
field direction reversal is equivalent to changing the sign 
of A. When the sign of A changes the potentials Vi 2(f) do 
not remain attractive for m > and so localized states of 
the Dirac excitations cannot be obtained. On the other 
hand for negative angular momentum states, m < — 1, 
Eq. (|14|) reveals that an attractive potential is possible 
when A < 0. Consequently, when the magnetic field di- 
rection is reversed one obtains localized states for the 
Dirac excitations for m < — 1. From the form of the ef- 
fective potentials in the Schrodinger equations of Xi an d 
X2, as given in the initial part of section m it can be eas- 
ily verified that if one changes the sign of A and replaces 
to by —(m + l) (for a negative to where to < —1) then xi 
and X2 just gets interchanged. Consequently we can say 
that the effect of the magnetic field reversal amounts to 
interchanging the components of the pseudospinors pre- 
viously obtained for the magnetic field pointing in the 
negative z-axis with their positive to values replaced by 
— (to + 1) where now to < — 1. It is interesting to point 
out that the interchange of the components of the pseu- 
dospinors can also be accompanied by some some relative 
phases between them. 

0.005 




up to a normalization constant and where y — Xr. The 
solution listed above is not square integrable and con- 
sequently it cannot be taken as the lower component of 
the pseudospinors. Consequently the lower component of 
the zero-energy states are all zero. The modulus squared 
pseudospinor components, V>i2 > for two specific cases 
as, to = 0, n — 2 and m = 1, n — 3, are plotted in Fig.0] 
and Fig.0 respectively. Here ■01*2 stands for either the jy 

upper component ip^ 1 '" 1 ^ or the lower component if;^'" 1 ^ 
of the pseudospinor for arbitrary values of the quantum 
numbers n and to. Both the curves show that the proba- 
bility densities of the excitations in both these cases are 
localized. 



FIG. 4: Plots of the normalized density profile of the upper 
component (solid) and lower component (dash) for n — 2, m = 
0. 



THE GROUND STATE SPECTRUM OF A 
MORE REALISTIC MAGNETIC FIELD 



Till now the present article discussed about a mag- 
netic field profile which becomes infinitely large at around 
r = 0. The advantage of choosing such a magnetic field 



6 




FIG. 5: Plots of the normalized density profile of the upper 
component (solid) and lower component (dash) for n — 3, m = 
1. 



is related to the fact that the Dirac-like excitations in 
graphene in the presence of such a field becomes localized 
and the spectrum of the excitations can be exactly solved. 
The exact solutions of the previous case could be found 
because one can match the problem of the graphene ex- 
citations in presence of an inhomogeneous magnetic field 
to the exact solutions of a 2D hydrogen atom bound 
state spectrum. In this section we will choose a magnetic 
field which produces a cut-off Coulomb scalar potential 
[Til [l5j | in the effective Schrodinger equation for the up- 
per component xi- I n this section we will mainly focus 
on the localized ground states of the Dirac excitations 
which can be exactly found. 

If one chooses the function f(r) as 



f(r) 



1 



a(r + a) 



a > . 



(33) 



where a is a positive constant whose actual dimension is 
that of length then the magnetic field is given by 



B.(r) = - 



1 



1 



a(r + a) (r + a) 2 



(34) 



In normal units the magnetic field will be obtained by 
multiplying the above expression by H/e. Clearly the 
magnetic field is a variable one and finite everywhere. In 
fact, it lies in the interval —2/ a 2 < B z (r) < 0. Let us 
now consider the sector m > 0. For the above choice of 
fir) the scalar potentials Vi^if) are given by 



Vi{r) = - 
V 2 (r) = - 



1 



2m + 3 
a(r + a) 
2m + 3 



a(r + a) (r + a) 2 



(35) 
(36) 



In the present case the two scalar potentials are differ- 
ent. It may be noted that apart from the constant 1/a 2 



which may be absorbed in the energy, V\{r) is the cut- 
off Coulomb potential albeit, in two dimensions. The 
other potential, V% (r) , is not a cut-off Coulomb potential 
since it contains an additional term. We would like to 
point out that the cut-off Coulomb potential of the form 
V(r) = — A/(r + a), where A > 0, has been extensively 
studied. The approximate solutions (l4l | as well as nu- 
merical solutions [151 ] of the Schrodinger equation in an 
cut-off Coulomb potential have been found using various 
methods. However in most of these studies the potential 
is treated in three dimensions where the angular momen- 
tum takes positive integer values starting from zero. On 
the other hand in the present case we are working in 
two dimensions and the coupling depends on m. Thus 
for m > — 1 the coupling is positive and the potential is 
an attractive one. On the other hand for m < — I the 
potential becomes repulsive. 

We note that Eq.© admits a e — solution and then 
it is given by 



Xx = r m+l/2( r . 



a)e 



-r/a 



(37) 



up to a normalization constant. Clearly \i i r ) — > as 
r — > oo, consequently Xi(r) is normalizable. 

The lower component of the pscudospinor solution of 
the Dirac equation can be written down once Eq.® is 
solved. The e = solution of Eq.© is given by 



X2 (r) = r-™- 1 / 2 



e r/a 

r + a 



(38) 



which is clearly non-normalizable for m > 0. Therefore 
the E — solution of the original problem is given by 



r m e ime (r + a)e~ r / a 




(39) 



where m = 0, 1, 2, • • •. The normalization constant is 



7V_ 



22m+3 



2-k a 2m+4 (m + 3) (2m + 3)r(2m + 2) 



(40) 



It may be noted that ip_ (r, 6) is normalizable for all 
values of the quantum number m > 0. This means that 
the e = state i.e, E = state is infinitely degenerate. 
Furthermore, there is no pseudospin down state corre- 
sponding to this level. For the case where m < 0, neither 
the solutions in Eq. ([3T]) nor Eq. (|3"5|) are normalizable. 
Consequently for m < all the states have non-zero en- 
ergy. The probability density of the zero-energy excita- 
tions for various angular momentum values are plotted 
in Fig. [Bl The curves show the localized nature of the 
states. 



A. Reversing the field direction 

For the sake of completeness one can also analyze 
the spectrum of the Dirac electron-like excitations of 
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graphene in the specific magnetic field as given in Eq. ([34]) 
with its sign reversed. In this case the vector potentials 
and the magnetic field are given by 

A x = -yf(r), A y = xf{r), (41) 

where /(r) is as given in Eq. (|33j) . The zero-energy solu- 
tions of Eq. ((6|) and Eq. ((9]) with the appropriate gauge 
potentials as given in Eq. (14"TT) are given as 

m+1/2 

Xi = — e r '\ X 2 = r-(™ +1 / 2 ' (r + a)e~ r / a , (42) 

r + a 

up to a normalization constant. From the above expres- 
sion it is seen that \i is n °t normalizable for any value 
of to. On the other hand, \2 is also not normalizable for 
to > 0. However, for to < the term — (to + ^) becomes 
positive and consequently \2 becomes normalizable for 

to = — 1, — 2, In other words when the direction of the 

magnetic field is reversed, the E — pseudospin down 
states are infinitely degenerate and there are no E = 
pseudospin up states. Explicitly the E = solutions are 
given by 

V^M) = N + y r -(m+l) e i(m+l)e^. + a y-r/a ) > 

(43) 

where m = — 1, —2, —3 • •• and N + is obtained from A_ if 
one replaces the m by — (m + 1). Except these replace- 
ment in A_ to obtain N + there can be a possible phase 
factor arising out of the magnetic field reversal about 
which we will discuss later. The important thing to note 
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FIG. 6: Plots of the normalized probability density of the 
zero-energy states vs y for different values of m — (long 
dot-dashed), m = 1 (dashed), m — 2 (solid) and m = 3 
(short dot-dashed). Here ip{ is the upper component of the 
pseudospinor for arbitrary values of m. 

is that when the magnetic field is reversed instead of the 
pseudospin up state we only have the pseudospin down 
state in the zero-energy solution. 



In this case also it can be noted that if one reverses the 
magnetic field and replaces to by — (m+1) (for a negative 
to where to < —1) then xi and xi just gets interchanged. 
The result of the field reversal essentially remains the 
same as discussed in the previous section. Magnetic field 
reversal amounts to interchanging the components of the 
pseudospinors previously obtained for the magnetic field 
pointing in the negative z-axis with their positive to val- 
ues replaced by — (to + 1). Presently to < — 1 and the 
interchange of the pseudospinor components can be ac- 
companied by a possible relative phase factor multiplying 
any of the components. 

We would like to note that many of these unusual spec- 
tral properties for 2D inverse-radial magnetic fields was 
also observed in the case of planar Schrodinger equation 
in the presence of a non-constant magnetic field [16j . In 
the present case the presence of pseudospin adds addi- 
tional terms to the effective potential and makes the sit- 
uation more interesting. 



V. DISCUSSION AND CONCLUSION 

The present article discusses about the possibility of 
localizing the massless Dirac excitations in graphene in 
the presence of some specific magnetic field profiles. The 
unique structure of the minimally coupled Dirac equa- 
tion for a massless particle in two dimensions gives an 
interesting physical and mathematical arena to look out 
for localized excitations. Localization of the excitations 
in graphene is an important topic as conventional elec- 
trostatic methods to localize excited modes fail. The 
inverse-radial magnetic field, which is directed along the 
z direction, has the potential to localize the graphene ex- 
citations near the origin. While there were many previ- 
ous works dealing with the topic of magnetic confinement 
of the excitations in graphene, most of those works used 
magnetic fields which vanished sharply at some radial dis- 
tance from the origin. The inverse-radial magnetic field 
defies this trend of choosing magnetic field profiles which 
can produce magnetic confinement of the excitations in 
the sense that the inverse-radial magnetic field does not 
sharply fall or rise at any distance from the origin. Al- 
though the magnetic field gradually falls away from the 
origin most of the excitations are localized near the ori- 
gin. 

It was pointed out in the article that in presence of 
the inverse-radial magnetic field the upper component of 
the pseudospinor of the graphene excitations follow an 
effective Schrodinger equation. An interesting observa- 
tion related to the effective Schrodinger equation of the 
upper component is related to the fact that this equa- 
tion turns out to be the Schrodinger equation of the 2D 
hydrogen atom. Consequently if one knows the bound 
state spectrum of the 2D Hydrogen atom one can also 
solve for the upper component of the graphene excita- 
tions in presence of an inverse-radial magnetic field. The 
lower component of the spinor can easily be found by the 
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application of the Dirac equation. Our work gives a nat- 
ural and physical application to the 2D hydrogen atom 
problem which by itself remains an esoteric and formal 
problem to deal with. 

In the case of the finite magnetic field the relevant 
Dirac equation for the excitations had a close connection 
with the 2D cut-off Coulomb potential. The Schrodinger 
equation for a cut-off Coulomb potential still remains 
to be analytically solved. Consequently for the realis- 
tic magnetic field the full excitation spectrum could not 
be presented because to present the full spectrum one has 
to solve an effective Schrodinger equation in presence of 
the 2D cut-off Coulomb potential. 

One of the interesting features of the spectra of the 
graphene excitations in presence of the inhomogeneous 
magnetic fields is related to the degeneracy pattern of 
the solutions. We have tried to point out some of the rel- 
evant degeneracy patterns and a general way how these 
degeneracies arise in the case of the inverse-radial mag- 
netic field. The ground state of the magnetically localized 
ground state turns out to be infinitely degenerate in both 
the cases. This is valid when the system size is infinite 
and the wave function vanishes at infinity. In principle, 
infinite number of electrons can be accommodated in the 
zero-energy states. For a realistic situation the degener- 
acy in the energy levels can not be infinity because of the 
finite size of the system. Fig. [3] shows that the position 
of the peak of the densities for various cases is shifting 
away from the origin as we increase the quantum num- 
bers (n,m). The distance of the point where the density 
maximizes can not be greater than the radius of the sys- 
tem. This implies that there is an upper cut-off in the 
(n, m) values. Therefore, the degeneracy is finite in a 
realistic situation. 

Throughout our analysis we have tried to see the ef- 
fects of a magnetic field direction reversal on the exci- 
tation spectrum. The result obtained shows up most 
prominently in the zero-energy states. When the mag- 
netic field is pointed in the negative z-direction only the 
upper component of the pseudospinor survives where as 
reversing the field direction kills the upper component 
and only the lower component of the pseudospinor ex- 
ists. In the most general case it was pointed out that 
a magnetic field reversal interchanges the pseudospinor 
components accompanied by a change of the m quan- 
tum numbers to —(m + f) where now m < — 1. It was 
also pointed out for both the cases, discussed in this arti- 
cle, such an interchange of the pseudospinor components 
could be accompanied by a relative phase arising out of 



the field reversal. In this section we can shed some light 
on the relative phase accompanying such a magnetic field 
reversal. 

It is interesting to connect these phenomenon related 
to the magnetic field reversal with the time-reversal like 
operation S as discussed in the review by Beenakker [2(| ■ 
Under the action of S the Hamiltonian of the excitation 
in Eq. (flj changes as 

SH(A)S~ 1 = H(-A) , (44) 
where A is the gauge potential giving rise to the mag- 
netic field. Although S is not the exact time-reversal 
operator, the exact time reversal operator changes the 
valleys, it acts like the time-reversal operator in a sin- 
gle valley. The action of S is like reversing the magnetic 
field direction. The operator S = iu^C where C is the 
complex conjugation operation. It can easily be verified 
that if one apply this operator on the pseudospinors one 
gets the pseudospinors of the reversed field case, as dis- 
cussed in this article. As a result of the S operation on 
the pseudospinors the component of the previous pseu- 
dospinor components gets interchanged with a relative 
sign change. This appearance of a sign can be easily ac- 
commodated in the scheme discussed previously where it 
was pointed out that the interchange of the pseudospinor 
components can be accompanied by a relative phase. 

To conclude this article it can be stated that the 
present article tries to analyze whether one can have 
magnetically confined excitations in graphene with in- 
homogeneous magnetic fields which do not have any fi- 
nite/infinite discontinuities in the non-zero finite range of 
the radial coordinate. The magnetic field profiles chosen 
turns out to be interesting because their specific func- 
tional forms allows one to use the results obtained for the 
bound states of 2D Schrodinger equations in the presence 
of a Coulomb or a cut-off Coulomb potential. It is ob- 
served that there exists an infinite number of localized 
excitations which are magnetically confined. In the fu- 
ture if one can produce a inverse-radial magnetic field 
then these exotic localized excitations can be observed 
experimentally. 
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